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Discovery of multi-phase cold accretion in a massive galaxy at z=0.7 
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ABSTRACT 

We present detailed photo+collisional ionization models and kinematic models of the multi- 
phase absorbing gas, detected within the HST/COS, HST/STIS, and Keck/HIRES spectra of 
the background quasar TON 153, at 104 kpc along the projected minor axis of a star-forming 
spiral galaxy (z = 0.6610). Complementary g'r'i'Ks photometry and stellar population mod- 
els indicate that the host galaxy is dominated by a ~ 4 Gyr stellar population with slightly 
greater than solar metallicity and has an estimated logM* = 11 and a logM v ; r = 13. Pho- 
toionization models of the low ionization absorption, (Mgi, Sin, Mgii and Cm) which trace 
the bulk of the hydrogen, constrain the multi-component gas to be cold (logT = 3.8 - 5.2) 
and metal poor (-1.68 < [X/H] < -1.64). A lagging halo model reproduces the low ioniza- 
tion absorption kinematics, suggesting gas coupled to the disk angular momentum, consistent 
with cold accretion mode material in simulations. The Civ and Ovi absorption is best mod- 
eled in a separate collisionally ionized metal-poor (-2.50 < [X/H] < -1.93) warm phase with 
logT = 5.3. Although their kinematics are consistent with a wind model, given the 2-2.5 dex 
difference between the galaxy stellar metallicity and the absorption metallicity indicates the 
gas cannot arise from galactic winds. We discuss and conclude that although the quasar sight- 
line passes along the galaxy minor axis at projected distance of 0.3 virial radii, well inside 
its virial shock radius, the combination of the relative kinematics, temperatures, and relative 
metallicities indicated that the multi-phase absorbing gas arises from cold accretion around 
this massive galaxy. Our results appear to contradict recent interpretations that absorption 
probing the projected minor axis of a galaxy is sampling winds. 

Key words: — galaxies: ISM, haloes — quasars: absorption lines. 



1 INTRODUCTION 

Over the last decade, simulations have shown that galaxy evo- 
lution is highly dependent on gas accretion occurring via two 
modes: hot and cold accretion. Current cosmological simulations 
demonstrate that the majority of gas accreted at early epochs 
onto galaxies occurs via the cold mode, which has temperatures 
of T ~ 10 4 - 10 s K and metallicities of Z < 0.01Z Q . Cold- 
mode gas is preferentially accreted along cosmic filaments/streams 
and have high densities and low cooling times providing a 
large suppl y of gas penetrating through hot halos surrounding 
galaxies (Keres et al. 2005; Dekel & Birnboim 2006; Ocvirk et al. 



20081: iKeres et all 120091: [Brooks et alj |2009t iDekel et alj 120091; 
Ceverino et al. 2010; Stewart et al. 201 la b; van de Voort et al. 
20111 : Ivan de Voort & Schavdl201ll ; lFaucher-Giguere et alj|20llfh 



It is expected that cold accretion should c omprise no more than 
7% o f the total Hi mass density at z ~ 1 l lKacprzak & Churchilll 

Unil)- 

It is further expected that cold accretion truncates when 
the host galaxy mass exceeds ~ 10 12 M Q , since infalling gas 
becomes shock heated to the halo viral temperature (~ 10 6 ) 
and is predicted to dramatic ally reduce the cold accretion cross- 
section to a tiny frac t ion (e . gjDekel & Birnboimf 2006 ; Kere s et all 
| 2005l: lOcvirk et al.1 120081; iDekel et alj 120091; IKeres et alj | 20~ ; 
Stewa rt et alj|201 lal: iBrooks et alj|2009l ; Ivan de Voort et alj|201 1 ; 



gkacprzak@astro.swin.edu.au 



van de Voort & Schayel 20111) of the observe d halo gas cross- 
section dKacprzak et alj2008l ; IChen et alj|20ld) . However, it is ex- 
pected that these dense filaments can still survive within hot halos 
and could provide an efficient means o f feeding massive galaxies 
with pristine gas (e.g. lKeres etai]|2005h . 

The study of absorbing foreground gas detected in back- 
ground quasar spectra allows us to probe these otherwise unob- 
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servable comic filaments and outflows. Mgn absorption is ideal 
for detecting cold mode and hot mode accretion, wind outflows, 
etc., since it probes gas with a large range of n eutral hydrogen 
colum n densities, 10 16 < N(Hi) <, W 1 cnT 2 dChurchill et al.l 
2000a]; iRigbv, Charlton. & Churchill 1 20021) with gas temperature 
of around 30,000-40,000 K and average total hydrogen den- 
sities of ~ 0.1 atoms cm dChurchill. Vogt. & Charlton] 1200 ll ; 



iDing. Charlton. & Churchilll l2005h . It has also been thoroughly 
demonstrated that Mgn absorption is produced within gaseous ha- 
los surrounding galaxi es and is not produced within the inter galac- 
tic medium (IGM) (see lChurchill. Kacprzak. & Steidell2005t) . 

Over the last decade, strong Mgn absorption has also been 
observed to directly trace 100-1000 km s~' galactic-scale outflows 
dTremonti et al.ll2007l; Iweiner et alJ|20Qg|: iMartin & Bouche1l2009l ; 
Rubin et al.ll2010l : ICoil et alJl201ll ; lMartin et alj|2012h that extend 
out to at least 50 kpc along the galaxy minor ax is dBordoloi et al.l 
1201 ll ; iBouche et alj|20lil ; IKacprzak et alj|2012l) . Galactic winds 
have been observed ov er a large range of r edshifts and detected us- 
ing a range of ions (see lSteidel et al.l2010l . and references therein). 
Correlations between galaxy colors and star formation rates with 
Mgn equivalent widt hs also indirectly suggest that absorption is 
produced in out flows dzibetti et alj|2007l : iNoterdaeme et al.ll20103 ; 
lNestoretal.ll201ll) . 

However, Mgn has been observed infalling (Mar tin et al.l 

120121) into highly inclined gala xies with velocities of 
100-200 km IT 1 dRubin et alj 1201 ll) . This is consistent with 
IKacprzak et al] J201 lbh who showed that absorption strength is 
correlated with the orientation of the galaxy major axis, implying 
that a significant fraction of weaker Mgn absorption systems are 
likely accreting toward the galaxy via cold flows. The bimodal 
azimuthal angle distribution of quasar sight-lines around Mgn ab- 
sorption selected galaxies also sugges ts that infall occurs along the 
proje cted galaxy major axis dBouche et al.l l20l"Tl ; IKacprzak et al.l 
|2012|) . These cold-flow streams likely produce a circumgalactic 
co-rotating gas component that is predominately infalling towards 
the galaxy and, in absorption, these structures are expected to have 
~ 100 km s velocity offsets relativ e to the host galaxy and in 
the same direction of galaxy rotation dStewart" etal. 2011b). These 
model s are consistent with p revious observations of ISteidel et al.l 
d2002h and IKacprzak et ail fcOlOal) that show Mgn absorption 
residing fully to one side of the galaxy systemic velocity and 
usually aligned with expected galaxy rotation direction, with the 
absorption essentially mimicking the extension of the galaxy 
rotation curve out into the halo. We expect low ionization states, 
such as Mgi, Mgn, Sin, Cn and Cm to be ideal for tracing cold 
mode accretion given metallicity, temperatures and densities 
expected. 

A reliable means of determining the origins of the absorbing 
gas is to obtain both the host-galaxy and absorption-line metallicity. 
Absorption-line metallicities for a handful of systems have been de- 
termined to range between [M/H] < -1.8 to -1 while existing near 



2004 ; Chen et alj 2005^ Tripp et al. 


2004 Cooksevetal. 


2008; 


Kacprzak et alj2010b; Ribaudo et al. 


l201lHThometalj2011 


). It is 



postulated that these extremely low metallicity absorption systems 
are likely accreting onto their host galaxies and possibly trace cold 
mode accretion, which is still expected for these sub-L* galaxies. 
In the rare case where absorption-line metallicities are larger than 
the host galaxy is suggestive that the absorption is probing winds 
jPerouxetal]|201ll) . 

Here we target a particular galaxy that has Mgn absorption 
consistent with disk-like kinematics, possibly tracing cold accre- 



tion. The absorption also contains a separate warm phase as in- 
dicated by separate strong Civ absorption that does not coincide 
with Mgn. We have obtained supplementary HST/COS data in or- 
der to determine the physical properties of the gas. In this paper, 
we perform kinematic and photo+collisional ionization models of 
multi-phase absorbing gas obtained from HST/COS, HST/ST1S, 
and Keck/HIRES, which is associated with star-forming spiral 
galaxy at z = 0.6610. In §|2] we describe our targeted galaxy and 
our data. We discuss our absorption-line analysis in §[3] In §[4] we 
describe the host galaxy properties determined from broad-band 
photometry and stellar population models. In §[5] we describe the 
results of our kinematics and photo+collisional ionization models 
and the physical properties of the absorbing gas. In § [6] we dis- 
cuss the possible origins of the absorption and our concluding re- 
marks are in §[7] Throughout we adopt an H = 70 km s~'Mpc~', 
Am = 0.3, Oa = 0.7 cosmology. 



2 TARGET FIELD AND OBSERVATIONS 

TON 153, also known as Q1317+227, is a bright (V=16.0 mag) 
quasar at z em =1.017. Inspection of a low resolution quasar spec- 
trum reveale d two Mgn absorption systems at z a bs = 0.29 and 
Zabs = 0.66 (Steidel & Sar genj|l992l) . Following a spectroscopic 
surve y of galaxies in clo se angular proximity to the quasar sight- 
line, ISteidel et alj j2002l) discovered galaxies Gl and G2 shown 
in Figure [TJ The quiescent early-type galaxy Gl has a redshift of 
Zgai = 0.6719 and the star-forming d isk galaxy G2 has a redshift of 
z g ai = 0.6610. lchurchill et all d2007l) demonstrated that Gl was as- 
sociated with a Lyo- complex that did not have any observable met- 
als even though it resides at D = 58.1 kpc; well within the 100 kpc 

I II 1 I — I 

where metals are expected (e.g., Chen et al. 2001; Kacprzak et al. 

2008; Chen et al. 2010; Tumlinson et al. 201 1). In a companion pa- 
per jChurc hill et al. 2 012bh . we further discuss Gl and its associ- 
ated absorption lines. The galaxy G2 is at D = 103.9 kpc and is as- 
sociated with extensive metal-line/LLS abso rption at z a ts - 0.6601 
dSteidel & Sargentlll992l;lBahcall et alll993lil99d;rChurchill et al] 
l2000al ; IDing, Charlton. & Churchilll 120051 ; Ichurchill et alj l2007h 



and is the focus of this paper. 



2.1 HST Imaging 

In Figure[TJwe present a 4700 second #Sr/WFPC-2 F702W image 
(PID 5984; PI Steidel) that was reduced using the WFPC-2 Associ- 
ations Science Products Pipeline 

(WASPlB- Note that the F702W 
filter provides a bandpass similar to a rest-frame Johnson B-band 
filter for galaxies at z ~ 0.6. Galaxy magni tudes and luminosities 
were obtained from lChurchill et alj d2012al) and are based on the 
AB system. The GIM2D software dSimard et alj |2002) was used to 
obtain q uantified galaxy morphologic al parameters that wer e pub- 
lished in lKacprzak et alj d2011bl) and lChurchill et alj d2012al) . 



2.2 Ground-based Imaging 

To further constrain the properties of G2, we analyzed multi-band 
g'r'i' and Ks imaging. The g'r'i' -bands were obtained using SPI- 
cam CCD imager on the APO 3.5 m telescope. SPIcam has a field 
of view of 4.78 x 4.78' and a spatial resolution of 0.14"/pixel. Our 



http://archive.stsci.edu/hst/wfpc2/pipeline.html 
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Figur e 1. — (left) Gr ound-based g'— (blue), i'- (green) and Ks- band (red) color composite centered on the quasar. Note that Gl is redder than G2 (see 
IChurchill et al]|2012bl) . Also note that the surrounding galaxies in the field are extremely red compared to both z ~ 0.6 galaxies, indicating that they are at 
much higher redshifts and not associated with the absorption, (middle) A 13" X 22" HST/WFPC-2 F702W image of the quasar field TON 153. The quasar is 
the brightest object in the northeast corner of the image. G2 is associated with a variety of metal-lines shown in Figure|2] A possible object seen within ~1" 
northeast of the quasar (G3) has not been successfully identified spectroscopically since the quasar is bright. Note the quasar sight-line passes near the minor 
axis of the moderately inclined disk of G2. — (top right) A Keck/LRIS spectrum of Gl containing Can H & K confirms the redshift of Gl to be z ga l = 0.6719. 
— (bottom left) A Keck/LRIS spectrum of G2 containing [On] emission confirms the redshift of G2 to be z ga l = 0.6610. 



observations were taken with on-chip binning of 2 x 2 which pro- 
vides a plate scale of 0.28"/pixel. The images were taken as part of 
a large survey program and were observed over 4 nights between 
March 2006 and March 2007 providing total exposure times in the 
g'r'i' of 5190, 4630 and 4350 seconds, respectively, with typical 
seeing of 1.1 - 1.6". 

Multiple frames were taken in each filter and each frame 
was individually reduced using standard IDL and IRAlfl pack- 
ages. Pixel-to-pixel variations were removed using a combination 
of dome and twilight sky flat fields. Due to PSF/seeing variations 
over the long exposures, cosmic rays were removed from each 
frame separately. The SPIcam pixels are sufficiently small that in- 
terpolation errors do not lead to significant photometric uncertain- 
ties. The astrometry was calibrated by matching field stars from 
each frame to USNO A2.0 catalog stars. 

The photometric zeropoints were established using a number 
of stars from the SDSS catalog. Since the APO g'r'i' filter/detector 
combination does not exactly match that of the SDSS survey, there 
can exist color terms between the two systems. We derived these 
color terms by analyzing the SDSS stellar photometry from over 
30 fields and produced a single photometric solution. We adopt the 
color term derived from the large set of observations and the aver- 
age corrections are of the order of 0. 1 magnitudes. 

Deep near-infrared (Ks) band images were obtained in 1994 

2 IRAF is written and supported by the IRAF programming group at the 
National Optical Astronomy Observatories (NOAO) in Tucson, Arizona. 
NOAO is operated by the Association of Universities for Research in As- 
tronomy (AURA), Inc. under cooperative agreement with the National Sci- 
ence Foundation. 



February 24 using the Kitt Peak 4-m Mayall telescope and the 
IRIM NICMOS 3 256 x 256 camera with 0.6"/pixel. These 
images were part of a more exten sive Mgn galaxy survey of 
ISteidel. Dickinson. & PerssorJ dl994l) . The images were reduced 
using the contributed IRAF package DIMSUMj. The photomet- 
ric zeropoints were establish ed using a number of stars from the 
2MASS point-source catalog dSkrutskie et al.ll2006l) . 

In Figure [T] we present a g'i'Ks color composite image cen- 
tered on TON 153 with galax ies Gl and G2 labeled. See companion 
paper lChurchill et alj (2012b) for further discussions of the proper- 
ties of Gl. Note that all other objects near the quasar sight-line are 
red compared to the two z ~ 0.65 galaxies, indicating that all the 
other objects are likely at much high er redshifts and t hat Gl and 
G2 are likely isolated objects. Churchill et al. (2012b) noted that 
the ROSAT X-ray luminosity limit for this field is four orders of 
magnitude less than what is expect for cluster centers and is con- 
sistent with the expected luminosity for early-type galaxy halos at 
z = 0.67. In addition, the absorption metallicities determined by 
IChurchill et al] J2012bl) for Gl, and those found here, are l-2dex 
lower than expected for inter-cluster medium. Thus, it is likely that 
the absorption does not arise in a cluster environment. 

The galaxy G3 is located within ~1" northeast of the quasar 
sight-line and has not been successfully identified spectroscopi- 
cally since the quasar is bright. It is possible that this object is 
responsible for either the absorption associated with Gl, G2, the 
z = 0.29 absorption system, absorption at the quasar redshift, or 



3 http://iraf.noao.edu/iraf/ftp/contrib/dimsumV3/. DIMSUM was con- 
tributed by R Eisenhardt, M. Dickinson, S. A. Stanford, & F. Valdez. 
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Figure 2. — Metal absorption-lines associated with the z = 0.6610 galaxy G2 obtained from HST/COS, HST/STIS and Keck/HIRES. In each panel the 
transition and instrument are labeled. The absorption-line data (black) are relative to the systemic velocity of G2. The regions highlighted by the tick-marks 
(red) indicate the detected absorption features and the location of the tick-marks are defined to be where the equivalent width per resolution element recovers 
to a 1 <r detection threshold. The sigma spectrum is shown below the data (green). 



othe r metal-lines identified in this sight-line. If G3 is at z = 0.67, 
than lChurchill et al. (2012b) estimates its mass to be a factor of ten 
less than the mass of G2 and would have a simila r metallicity to 
G2 ac cording to mass-gas metallicity relations (e.g.. lSavaglio et all 
2005). Thus, if G3 is at the same redshift of G2, it could be consid- 



ered as a satellite of G2 and would occupy the same gaseous and 
dark matter halo. 

Photometry for calibration and s cience was extracted using 
SExtractor teertin & Arnoutsl [l996l) using the MAGAUTO mea- 
surements. Corrections for Galactic dust extinction were a pplied to 
the galaxies using the dust maps of Schle geTet alj dl998l) . We ob- 
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tained final dust-, color- and seeing-corrected AB magnitudes for 
Gl of m g > = 23.43 ± 0.03 = 22.17 ± 0.03, wy = 21.01 + 0.03 
and m Ks , = 19.4 ± 0.1 ^Churchill et"al]|2012bh and for G2 of 
m s , = 23.23 ± 0.03, m r , = 22.32 ± 0.03, m t = 21.41 ± 0.03 and 
m Ks , = 19.9 + 0.1. 



2.3 Galaxy Spectroscopy 

The galaxi es shown i n Figu re [T] were spectroscopically iden- 
tified by ISteidel et al. 1 2003) a nd their spectra were first pre- 
sented in IChurchill et alj J20071) . A possible object seen within 
~1" from the quasar seen in the image has not been success- 
fully identified spectroscopically since the quasar is bright. The 
de tails of the Keck/LR IS s pectroscopic observations can be found 
in 



tails ot the Keck/LRls> spectroscopic observations can be round 
Steid el et alj fe002h and lChurchill et al.1 ( 120071) ; IChurchill et al.1 



j2012bh . The spectra are both vacuum and heliocentric veloc- 
ity corrected. Galaxy Gl, identified by Can H & K absorption 
yields a z ga i = 0.6719, is associated with broad Lya- complex 
that s pans 1400 km s -1 and yet contains only v ery weak metals 
lines ( Churchill et a lj200l ; IChurchill et al.l 2012b). Galaxy G2 was 
identified by a [On] emission line, placing it at z gl ,i = 0.6610. This 
galaxy is associ ated with extensive metal - line/LLS absorption at 
z abs = 0.6601 teahcalletal.lll993lll996l:IChurchilletalJl2000j ; 



iDing. Charlton. & Churchillll2005l ; IChurchill et alj2007h and is the 
focus of this paper. We discuss G2's associated metals-lines in the 
next sections. 



2.4 Quasar Spectroscopy 

Some of the z„bs = 0.6601 absorption properties were measured 
from a 3600 second Keck/HIRES (R ~ 45, 000) exposure observed 
in 1995 January. The details of the observ ation and data reduction 
are described in IChurchill & Vogtl ( 1200 ll) . In addition, a 12,000 
second HST/STIS E230M (R ~ 30, 000) exposure was obtained 
(P IP 8672; PI Churchill) and the dat a reduction details are found 
in lDing. Charlton. & Churchilll J2005h . 

Recent HST/COS observations of the quasar (R ~ 18,000) 
were obtained using the FUV G160M grating and the NUV G185M 
grating (PID 11667; PI Churchill). The FUV observations were 
centered at 1600 A and took place on 25 June 2010 with an ex- 
posure time of 12,580 seconds. The NUV observations consisted 
of two exposures that occurred on 26 May 2010 and were opti- 
mally co-added. The first 5420 second NUV exposure was cen- 
tered at 1921 A and the 4970 second NUV exposure was centered 
at 1941 A. The COS spectra were reduced using the standard HST 
IRAF pipeline. All spectra are both vacuum and heliocentric veloc- 
ity corrected. 

It is important to note that the wavelength solutions across 
these multiple instruments are consistent at the sub-pixel level. This 
was verified across all spectrographs by centroiding common ion- 
izati on absorption lines from the data presented here and the data 
from lChurchill et alj d2012bh . 

Analysis of the absorption profiles was performed using inter- 
active software (see lChurchill et al.ll999ll2000al ; IChurchill & Vogtl 
l200ll) for local continuum fitting, objective feature identification, 
and measuring absorption properties. Velocity widths and equiv- 
alent widths of the absorption systems are measured between the 
pixels where the equivalent w idth per resolution e lement recovers 
to the 1 cr detection threshold dChurchillet alll999l) . 
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Figure 3. — (top-left) COS/HST Ovi il032 transition ass ociated with G2 
is blended with Ly/ ? from the z =0.67 Gl Lya complex (Churc hill et alj 
120071 ; IChurchill et alj|2012bl) . The Ly/3 fit (red) was constrained by simul- 
taneous Voigt profile fits to the Lya, Lyy, Lye, and Lyf. — (top-right) 
COS/HST Cm /1977 transition assoc iated with G2 is blended with a Lyy 
from t he z =0.67 Gl Lya complex JChurchill et"ai]|2007l ; IChurchill et al] 
l2012bl) . The Voigt profile fit was constrained by the Lya, Ly/3, Lye, and 
Lyf which produce the Lyy fist shown (red). — (bottom-left) The de- 
blended Ovil032 line associated with G2. — (bottom-right) The deblended 
C in /1977 line associated with G2. 



3 QUASAR ABSORPTION-LINE ANALYSIS 

In Figure [2] we present the absorption-line data and in Table [JJ 
we present the measured equivalent widths and column densi- 
ties. We discuss below how the column densities and/or limits 
were computed for each transition. We further discuss how we ac- 
count for higher redshift hydrogen-line blends that contaminate the 

Ovi /11032 and Cm A977 transitions (see Figure[3]l. 

The Hi column density was adopted from Churchill et al. 

( 120071) . obtained by simultaneously fitting the Lya, Ly/?, and Ly- 
man break obtained from a FOS spectrum, was determined to be 
log[W(Hi) ] = 1 8.3 ± 0.3. This is consistent with the results of 
iRao et all |2006), who fit only the Lya line, and determined the 
column density to be 18.57 ± 0.02. If we were to use the Rao et al 
value, our derived metallicities would decrease be 0.27 dex. 

We also note upon analyzing the Civ /1/11548, 1550 doublet, 
we have deemed the reddest component seen in the C iv A 1 55 1 tran- 
sition, at roughly 80 km s -1 , to be real and has been included in our 
analysis. This feature does not appear in the Civ /(1548 line since 
this region is contaminated by poor s ky subtraction. This compo- 
nent w as omitted by the analysis of IDing. Charlton. & Churchilll 
(2005) but included here. 



3.1 Accounting for Blends in O vi ,11032 and C m ,1977 

IChurchill et"al] d2012bh performed detailed Voigt profile (VP) fits 
to the Lyman series, using the HST/COS spectra, in order to de- 
termine the gas-phase properties of the Gl absorption complex. 
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Table 1. Absorption Properties" 



Ion 



Instrument/ 
Telescope 



W r 
(A) 



\ogN(X) 



W blm 
(A)* 



h \ogN(Xf 



W r ed b {ogN(X) n 
(A)* 



Mgn ,12796 
Mgn ,12803 
Mgll 

Mg 1,12853 
Mgl 

Sii ,12515 
Sii ,11846 
Sii ,11631 
Sii ,11562 
Sil 

Sin ,11527 
Sill 

Si iv A 1394 
Siiv ,11403 
SUV 

Cn ,11036 
CH 

Cm ,1977 

cm 

Civ ,11548 
Civ ,11551 
CIV 

Ovi ,11032 
Ovi ,11038 
OVI 

HI 



HIRES/Keck 
HIRES/Keck 



0.348±0.007 
0.198±0.006 



HIRES/Keck 0.039±0.005 



HIRES/Keck 
STIS/HST 
STIS/HST 
STIS/HST 

STIS/HST 

STIS/HST 
STIS/HST 

COS/HST 

COS/HST 

STIS/HST 
STIS/HST 

COS/HST 
COS/HST 

FOS/HST 
COS/HST 



<0.0066 
<0.053 
<0.019 
<0.023 

0.067±0.013 

<0.019 
<0.023 

>0.018±0.005 c 

0.467±0.008 

>0.508±0.027 rf 
0.417±0.026 

0.233±0.012 
0.105±0.010 



13.11+ 



11.54+ 



<11.8 



13.16+ 



<12.4 



13 39+0.09 



14 20+° 03 



14 41+0 .06 
1H - -0.05 



1 4 49+" " } 
1H - -0.03 

18.3±0.3 C 



0.348±0.007 
0.198±0.006 

0.039±O.OO5 

<0.0066 
<0.053 
<0.019 
<0.023 

0.067±0.013 

<0.019 
<0.023 

>0.018±0.005 c 

0.431 ±0.007 

0.315±0.024 
0.189±0.021 

0.189±0.011 
0.076±0.008 



13 H+ uu/ 
-0.07 

11 54+0.05 
11 ■-'^-0.05 

sSll.8 
13 16+ 011 

<12.4 



> 13 39+° 09 



14.18*$ 



14 11+° 07 
1H - -0.05 



14 33+o «» 
18.3±0.3 e 



<0.0049 
<0.()060 

<0.0049 

<0.0066 
<0.053 
<0.019 
<0.023 

<0.018 

<0.019 
<0.023 

<0.012 

0.036±0.004 

>0.193±0.031 rf 
0.228±0.018 

0.044±0.006 
().029±0.007 



sSll.06 
<10.58 

^1 1.8 
<12.6 

^12.4 
<13.1 



12.88*$ 



14 15+ 008 



13.94+ 



16 03+° 18 



" The equivalent width limits are computed at the 3<x-level and are for an unresolved cloud. The column densities are measured using the AOD method and 
the column density limits are computed using the COG analysis (see text). 

* W r and logN(X) are quoted for absorption blueward (W bl " e ) redward (W r ed ) of the galaxy systemic velocity. The velocity windows used to compute 
logA'(X) are —240 < v a b s < -3 km s for gas blueward of the G2 systemic velocity and —3 < i>„fts °S 220 km s for gas redward of the G2 systemic velocity. 
c Cn is blended with an unknown line at v a / JS < — 100 km s - ' , which is why the Cii is expressed as a lower limit. The exact value quoted is valid over the 
velocity window shown in Figure|2] 

d Crv il551 is affected by poor sky subtr action and is expressed as a lower limit.The exact value quoted is valid over the velocity window shown in Figure|2] 
e The Hi column density is adopted from lchurchill et all |2007[) and was determined by simultaneously fitting the Ly a, Ly/3, and the Lyman limit absorption. 



The VP fits were performed using our own software M INFIT 
dChurchill & Vogtl 1200 ll ; IChurchill. Vogt. & Charltor] 1200 ll) . The 
HST/COS instrumental line spread functions for the spectrograph 
settings, and also for the observed wavelength of each modeled 
transition, was computed by int erpolating the on-line tabulated data 
jDixon et alj2O10l ; lKrissll2Ol ll) . 

Churchill et al. (2012b) determined that a portion of the Gl 
Ly/3 complex was blended with the Ovi ,11032 associated with G2, 
while a different portion of the Gl Lyy complex was blended with 
C 111,1977 associated with G2 as shown in Figuref3](the full Ly/3 and 
Lyy complex is not shown since it spans roughly 1400 km s" 1 ). 

The VP fits of the Ovi ,11032 blended portion of the Ly/3 com- 
plex was well constrained by the non-blended Lyo-, Lyy, Lye, and 
Ly ( lines. The VP fits of the Cm ,1977 blended portion of the Lyy 
complex was well constrained by the non-blended Lya, Ly/3, Lye, 
and Lyf lines. We show the blended and deblended Ovi ,11032 and 
Cm A911 lines in Figure [3] Given that Gls hydrogen lines are un- 
saturated and can be well modeled using most of the Hi series, we 
are confident that we have a robust correction for the spectral shape 



of the Ovi ,11032 and Cm ,1977 lines. Furthermore, these correc- 
tions only effect regions redward of the galaxy systemic velocity. 



3.2 Apparent Optical Depth Method 

We employed the apparent optical depth method (AOD) to mea- 
su re the column densi t ies for eac h transition using t he for malism 
of ISavage & SembacH Jl99ll) and lChurchill & Vogtl feOOll) . In the 
cases that multiple transitions of a given ion were measured, we 
computed the optimal weighted mean column density in each ve- 
locity bin. Since the weighted mean requires the inverse square of 
the uncertainties, in these cases we treated the non-normal distribu- 
tions in the uncertaint i es of t he optical depths using the quadratic 
model o flD'Agonstinil J200C|) ; |D' Agonstini & Rasol?200rj|) : lBarIowl 
d2003l) . The quadratic model approximates the probability distribu- 
tion of the asymmetric uncertainties in each optical depth data point 
by a parabola fit to the dimidated Gaussian constructed from the 
upward and downward measured optical depth uncertainties. Rela- 
tive to a treatment that neglects the non-normal distribution in the 
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Figure 4. — The doublet optimally combined apparent optical depth distri- 
bution for Ovi (orange), Civ (black) and Mgn (red). The Hi (blue) optical 
depth distribution was produced by combining the entire Lyman series. The 
total column densities are integrated over the velocity bins. Note that the 
Hi is still saturated in the center of the profile. Also note the velocity struc- 
ture difference between the profiles and that the majority of the hydrogen is 
at velocities consistent with the Mgn whereas the majority of the Civ has 
velocities consistent with a small fraction of the total Hi. 
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Figure 5. — Color-colo r diagram for g - r versus — K. Lines are 
iBruzual & Charloj 12003^ stellar population models for an exponential star 
formation history and a|Chabrier ( 2003) initial mass function. Stellar metal- 
licities and formation redshifts are also labeled. The points (magenta) are 
galaxies Gl and G2, where G2 is the fo cus of this pape r and d iscussions of 
Gl can be found in a companion paper Churchill et al. 1 2012b). A redden- 
ing vector for E(B - V) = 0. 1 at the rest-frame of the galaxies is provided 
for reference. G2 is consistent with stellar population model of solar metal- 
licity. 



uncertainties, the method effectively provides a point-by-point bias 
correction to the resulting mean. 

In Figure|4]we present the optimal column density distribution 
for each transition. The Mgn, Civ, and O vi AOD distributions were 
computed using the doublet. For the Civ A 1548 line, we masked 
out pixels that were unusable as indicated by the two spikes in the 
sigma spectrum visible in Figure [2] Note the difference in the op- 
tical depth as a function of velocity of each transition, especially 
between Mgn and Civ. We also employed the AOD method to com- 
pute the remaining column densities presented in Table [T] 

Even though the Lyman series is saturated, we have optimally 
combined all the series lines to produce the Hi column density dis- 
tribution presented in Figure [4] Note that some central pixels are 
saturated and we therefore can not obtain a total column density, 
however, we can compute the Hi column density over a range of 
velocity windows where there is no saturation. 



3.3 Curve of Growth Analysis 

We use the curve of growth analysis (COG) to determine the col- 
umn density limits of various transitions of silicon and for Mgn, 
Mgi and Cn redward of the galaxy systemic velocity. The equiva- 
lent width limits quoted in Table [T] are 3cr limits for a unresolved 
single cloud. The single cloud assumption for silicon is motivated 
by the Sin /11527 absorption-line, given that we would likely expect 



to detect the remaining silicon transition as a single cloud compo- 
nent consistent with the strongest Mgn cloud located at roughly 
-180 km s- 1 . 

The measured equivalent width limits are small and being 
on the linear part of the curve of growth implies that the column 
densities are mostly independent of the Doppler parameter (/?). In 
fact, for b > 5, the predicted column density are independent of 
b. We find that the Sin column density is best constrained by the 
Sin /12515 data and the Siiv column density is best constrained by 
the Siiv ,11394 data. 

There also exists HST/FOS spectra of this quasar that cov- 
ers additional silicon transitions, among other metal-lines, however, 
we find that the equivalent widths were inconsistent with the pre- 
dictions by the COG analysis; this was not the case for the STIS 
or COS data. We n o te that the FOS measu rements published in 
Jannuzi et aD dl998h : IChurchill et all f2000a) are plagued by unre- 
solved blends from the other Hi lines and therefore we do not use 
the FOS data in our analysis. 

It is important to state that although the curve of growth anal- 
ysis seems to syst ematically under estimate the column densities 
of the metal-lines (Prochaska 2006), the silicon transitions do not 
provide sufficient constraints on our analysis and do not affect our 
results. We provide these column density measurements for com- 
pleteness. 
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Figure 6. — (top) Overlay of the Lya, Mgii A2796, Cm ,1977, Civ il551 
and Ovi ^ 1032 transitions. Note that the low ionization Mgii /12796, 
Cm /1977 trace the bulk of the hydrogen while the higher ionization states 
trace both low and high hydrogen column densities. The line morpholo- 
gies of Civ and Ovi are quite different to those of Mgii and Cm. — (mid- 
dle) AOD column densities of Mgii (blue circles), Cm (black squares), Civ 
(green triangles) and O vi (red pentagons) as a function of different velocity 
bins. The data points are located in the center of the bin, while the shad- 
ing indicates the full velocity width of the bin. — (bottom) Same as above 
except the values are normalized to JV(Civ). Note that Civ and Ovi have 
similar behavior blueward of the systemic velocity. 



4 RESULTS: GALAXY PROPERTIES 

In Figure Q] we present a g'i'Ks color composite centered on the 
quasar with galaxies Gl and G2 labeled. We find that G2 tends to be 
bluer than Gl, which likely indicates a younger stellar population. 
This is coincident with the fact that Gl exhibits weak (3 A) [On] 
emission. 

In Figure [5] the g - r versus i — K c olors of Gl and G2 are 
compared to the lBruzu al & Chariot 1 2003 ) stellar pop ulation mod- 
els obtained using EzGal ( Mancone & Gonzalezl2012l) for an expo- 
n ential star forma tion history, with an e-folding time of 1 Gyr, and 
a lChabried ( l2003l) initial mass function. The K-band is crucial in 
breaking degeneracies between other star formation history models 
since it is very sensitive to the old stellar populations. We tested a 
variety of models and find that anything but a prolonged r = 1 Gyr 
exponential model do not fit t he data. 

In our companion paper IChurchill et al. (2012b) we discuss 
Gl, although we note that Gl is redder, formed at an earlier epoch 



and is more massive. The model comparison suggests that G2 is 
dominated by a ~ 4 Gyr stellar population with slightly greater 
than solar metallicity abundance and formed at redshift z ~ 2. We 
also note while G2 has an [On] rest equivalen t width of 3.0 A, it i s 
consistent with galaxies having U - B > 1.0 ^Cooper et alj|200of) : 
red galaxies with little-to-no star formation. Thus, the stellar pop- 
ulation of G2 should be dominated by an older population as the 
model suggests. Adopting the A"-band mass-to-light ratio associ- 
ated with this stellar population model, we estimate the total stellar 
mass of G2 to be M, ~ 1 x 10" M G . 

We further estimate the halo virial mass from the stellar mass 
using GalMass ( Stewart 2011). GalMass uses abundance matching 
models from lMoster et al.l " 2010^ along with semi-empirical fits to 



observed galaxy g as fractions to conve rt between the stellar and gas 
mass described in JStewart e t al. 2009). There is roughly a 0.25 dex 
uncert ainty in M ylI at a fixed M, due to the s ystematics in estimates 
of M, teehroozi. Conrov. & Wech sler 2010) . The mean ga l axy gas 
and baryonic masses are also estimated from lStewart et al] ( 201 la) 
who employed the baryonic Tully-Fisher relat ion jMcGaughl2005h 
the stellar, gas, and dynamical mass relation (Erb et al. 2006), and 
galaxy-gas fraction and stellar mass relation dStewart et al.l 2009). 
We estimate the halo virial mass, gas mass and baryonic mass to be 
log M vlr = 12.9, log M m = 10.1 log M b = 11.1, respectfully. 



5 RESULTS: ABSORPTION 

From Figure[2]and Figure[4] it is clear that the absorption properties 
vary significantly as a function of ionization level. In Figure[6] we 
overlay selected transitions to further demonstrate this. It is clear 
that Mgii and Cm (along with Sin and Mgi) have similar kinemat- 
ics and trace the bulk of the hydrogen. While Civ and Ovi trace 
some of the same gas, their kinematics and their relative absorp- 
tion strengths differ. Furthermore, they both exist where there are 
diffuse hydrogen column densities and no Mgii absorption. 

In the bottom panels of Figure [6] we show the AOD column 
densities as a function of the velocity for each transition. Note that 
the Ovi follows similar abundances to Civ across the profile, while 
Mgii and Cm differ and decrease toward the positive velocities. It 
appears quite clear that the kinematics and the abundance patterns 
show that Civ and Ovi trace different phases of gas and this differ- 
ence become greater redward of the galaxy systemic velocity where 
the hydrogen column density decreases by 2 dex. 

Given the observed kinematic and abundance profile differ- 
ences, we have broken the absorption regions into two separate 
components and model them independently. In Table [T] we list 
the total column densities for gas blueward, N(X) bl " e , and redward, 
N(Xy d , of the galaxy systemic velocity. In the following subsec- 
tions, we present models to explain the observed absorption kine- 
matics and also apply photo+collisional ionization models to deter- 
mine the origins and the physical properties of the gas. 



5.1 Kinematic Models 



ISteidel et al] J2002h obtained the rotation curve for G2 and is rep- 
resented in Figure [7] Note that all of the Mgii, and the majority 
of the Cm, resides to on e side of the galaxy systemic velocity. 
IChurchill & Steidell l2003h noted that the C iv, and here we see the 

Ovi, resides on both sides of the galaxy systemic velo city. 

T o explore the origins of the Mgii absorption, Steidel et al. 
d2002l) employed a simple lagging halo model that extended the 
galaxy rotation velocity out into a co-rotating gaseous disk/halo. 
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Figure 7. — (left) The data points show the rotation curve of G2 obtained from lSteidel et alj J2002l) . Below are a selection of absorption transitions where the 
hashed region (blue) shows the rotating thick disk modeled velocities, while the solid shading (red) indicates the wind model velocity predictions. Note that 
all of the Mgn, and the majority of the Cm, resides to o ne side of the galaxy systemic velocity while Civ, and some Ovi and Cm also resides redward of the 
galaxy systemic velocity. — (right) The wind model of Bouche e t al.l 1201 ll) applied to G2. — (top-left panel) The wind modeled cone, having an opening 
angle of 30 degrees, viewed face-on where the galaxy is represented as the filled ellipse and with the QSO line-of-sight marked. — (top-right panel) A side 
view. — (bottom-left panel) The average cloud line-of-sight velocities as a function of position within the wind. The QSO location is represented as the filled 
circle. — (bottom-right panel) The line-of-sight velocity distribution of the clouds at the location of the quasar. We find that an outflow velocity of 150 km s _1 
produces a good match to the data and comparable to the results found by Bouche et alj l20l"ll) . Note the peak of the density distribution coincides with the 
bulk of the Civ absorption. 



They determined that a lagging halo model can account all of Mgii 
absorbing gas kinematics, even though the absorption occurs at 
D = 103.9 kpc. The predicted model velocities are highlighted by 
the hashed (blue) region below the galaxy rotation curve in Fig- 
ure [7] Further note that the observed rotation curve does not go 
deep enough to reach the flat part of the curve. This implies there 
is additional/higher galaxy kinematics that were not included in the 
model, thus the vertical solid (blue) line on the left of the absorption 
profiles would move towards bluer velocities. It is quite clear that 
the low ionization gas seems coupled t o the disk kinematics simi- 
lar to the models of IStewart et alj J201 lbl) . However, the observed 
co-rotating Men abso rption extends to larger D than the models of 
IStewart et alj ( feOllbl) since, at z = 1 .4, their simulated disks extend 
only to ~40 kpc. Thus, it remains unclear if there is size evolution 
as a function of redshift or if the simulations (without metals) prop- 
erly trace the metal-lines or if this is a result of simulating only a 
few galaxies. Aside from these caveats, the similarities between the 
observations and the simulations are suggestive that this absorption 
system exhibits signatures of cold mode accretion. 

Note that the rotating disk model does not account for the all 
of the absorption. Perhaps that given that the quasar line-of-sight 
passes along the minor axis of the galaxy and that G2 is undergo- 
ing some star-formation (indicated by the [On] emission shown in 
Figure[TJi, then it is possible that some of the absorption arises from 
outflows. 

By analyzing 10 z~0.1 Mgii absorbers selected from 



ItCacprzak et alj d201 lah . lBouche et alj J20 1 lh found a bi-modal dis- 
tribution of the azimuthal orientation of the quasar sight-lines: half 
of the sight-lines aligned with the major axis and half within a =30 
degree of the minor axis. The bi-modal azimutha l angle distribu- 
tion was later confirmed bv lKacprzak et alj d2012h using a sample 
of 88 absorption-selected galaxies (W r (2796) > 0.1 A) and 35 non- 
absorbing galaxies (W r (2796) < 0.1 A). These results indicate that 
both gaseous disks and strong bipolar ou tflows could contribute to 
Mgii cross-section. iBoucheetaflbOllh also applied a simple bi- 
conical wind model that was able to reproduce the observed Mgii 
kinematics for the sight-lines aligned with the minor axis. We ap- 
ply their model here in an effort to reproduce the observed absorp- 
tion kinematics. Their galactic wind model consists of 10 5 "clouds" 
moving at a constant velocity (V m „) that is determined from the ab- 
sorption data. The clouds are contained within a cone that has a 30 
degree opening angle. 

Figure [7] sh ows the wind model fo r G2 using the orientation 
parameters from Kacprzak et al.1 i20l lbl ). The top-left panel shows 
the cone view face-on and the top-right panel shows a side view of 
the cone, where the galaxy is represented as the filled ellipse and 
with the QSO line-of-sight marked. The bottom-left panel shows 
the average cloud line-of-sight velocities a function of position. The 
bottom-right panel shows the distribution of the cloud line-of-sight 
velocities along the quasar sightline. The wind speed is tuned to 
match the observed velocity range. For this particular case, we find 
that an outflow velocity of 150 km s _I produces a good match to 
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the da ta and are comparable to the results found bv lBouche et al.l 

fconh . 

In the lower-right panel, we see the predicted absorption dis- 
tribution peaks at roughly 50 km s~ x , where the majority of the Civ 
and some Ovi resides. The wind model velocity range is also shown 
as the solid shaded region over the absorption profiles. The model 
could also account for the observed Civ and Ovi blueward of the 
galaxy systemic velocity, although the model predicts only predicts 
of few percent of the gas is expected at these velocities. 

If the opening angle or the wind speed is increased, the model 
velocity range would also increase and could include all of the ob- 
served absorption. However, increasing these parameters will not 
reproduce the observed optical depth distribution of the cold and 
hot gas. While our current model is tuned to reproduce the absorp- 
tion reward of the galaxy systemic, a different model would not 
be able to predict the strong absorption residing at -200 km s . A 
model with increased opening angle and/or wind speed would also 
predict additional absorption beyond +200 kms~' where none is 
observed. Therefore, it is very difficult to reproduce both the ve- 
locity and optical depth distributions of the cold and hot gas using 
only the wind model. 

The wind model can not account for the bulk of the Mgn, Mgi, 
Sin and Cm. Thus, although the galaxy is forming stars, winds are 
likely not responsible for the cool gas, but could be responsible for 
the hot gas. One would naively expect that the infalling gas would 
be metal poor while the outflowing gas to be metal enriched. We 
explore the gas-phase metallicities in the next subsections in order 
to determine the possible origins of the absorption. 

5.2 Cold Gas Phase 

To determine the physical properties of the c ool gas blueward of 
the galaxy systemic velocity, we use Cloudy [Ferlan dl-Talll 19981) 
to model the ionization conditions. We follow the standard assump- 
tion of a photoionized uniform sl ab of gas that is in ioniza tion equi- 
librium and is illuminated with a Haardt & Madau] J20 1 lh ionizing 
spectrum, where the UV photons arise from quasars and galaxies. 
The ionization parameters, U, and the metallicity of the gas are 
varied to match the observations of N(X) bl " e in Table[T](column 6). 

In Figure [8] we show the model results computed for 
log[iV(Hi)] = 18.3 ± 0.3 and a metallicity [X/H] = -1.67. We 
find a narrow range of ionization parameters that reproduce the 
cool gas phase. The thin curves show the models while the thicker 
curves show the model values permitted by the data. The Cm 
provides the tightest constraints on the ionization parameter of 
-3.25 < log(7 ^ -3.21. Less stringent constraints are also found 
from Sin and Mgn yields -3.38 < logf/ s: -3.16. 

In Figure [8] the same Cloudy models are shown as a function 
of U and metallicity. The constraints placed on U and the models 
confine the allowed range of metallicity. The Cm measurements 
limit the metallicity to -1.68 < [X/H] «S -1.64. Recall that the 
Cm column density measurement blueward of the galaxy systemic 
velocity is the unblended portion of this transition. If we were to 
include the small contribution of the Cm redward of the galaxy 
systemic velocity our results would still be consistent. The Sin and 
Mgn yields a less stringent of -1.73 < [X/H] < -1.63. Thus, 
this cool gas component has low metallicity and a low ionization 
parameter. 

In summary, the Cm column density provides the tightest con- 
straint and yields a log(/=-3.23 ± 0.2 and [X/H] = -1.66 ± 
0.02. The gas is primarily ionized since the ionization fraction 
X(Hi)=jV(Hi)/W(H)=0.088 + 0.002 with a hydrogen number den- 
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Figure 8. — (top) The Cloudy predicted column densities, as a function 
of the ionization parameter, for a metallicity of [X/H] = -1.67 and a 
log[A'(Hi)]=18.3. The bold portions of the curves show the observed col- 
umn densities, N(X) bh " blueward of the galaxy systemic velocity (Table[T| 
column 6). The vertical lines show the ionization parameter constraints pro- 
vided by the data. The Cm provides the tightest constraints of -3.25 < 
logj/ < —3.21 indicated between the solid vertical lines (gold region). Sin 
and Mgn provide secondary constraints of —3.38 < logU < -3.16 as noted 
by the vertical dashed lines (yellow region). — (bottom) The ionization pa- 
rameter as a function of metallicity for the above model. The shaded regions 
show the model ranges due to the measurement errors in N(X) bl " € . The ver- 
tical lines are the ionization parameter constraints provided by the top panel 
and provides the region of acceptable metallicities. The overlapping shaded 
regions of each ion within the constrains provided by U have a solid shad- 
ing and outlines the allowed metallicities of this cold gas. Again, Cm places 
the tightest constraints on the metallicity —1.68 < [X/H] <, —1.64. Thus, 
the cool gas has low metallicity and a low ionization parameter. 
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Figure 9. — (top) The photo+collisional ionization models 
I Churchill & Klimek 20121) as a function of temperature and hydro- 
gen density. The lines indicate abundance ratios which are independent 
of the hydrogen column density. The solid line is the Civ/Ovi column 
density ratio along with its measured error (dotted lines). The long and 
short dashed lines are the Chi/Civ and Cm/Ovi column density ratio, 
respectively. The Cm/Civ and Cm/Ovi dashed curves are lower limits on 
the allowed temperature range for a given hydrogen density. Note that full 
self consistency is when the dashed curves to fall below the black curve. 
The gas has a hydrogen density of log;i H > -3.5 and the Civ/Ovi provides 
a temperature constraint of log? = 5.23 - 5.29. — (bottom) The metallicity 
as a function of hydrogen density for various hydrogen column densities. 
We provide conservative limit of logA'(Hi) < 16.8 [see text] which yields 
[X/H] < -0.5. 



sity of log(n H )=— 1.85 ± 0.02. The physical size of the cloud is 
L = N H /n H =521 ±2 pc. 

Note that Civ and Ovi do not appear on these plots. The col- 
umn densities predict for this photo-ionized modeled gas would be 
logiV(Crv) = 11.88 - 11.94 and logJV(Ovi) = 7.29 - 7.35. This 
predicted column densities are no where near the measured values 
of logN(Civ)=14.11 and logN(Ovi)=14.33. Thus, it is likely that 
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Figure 10. — The phot o+collisional ionization models obtained from 
IChurchi ll & Klimek] \20\% as a function of temperature and column den- 
sity normalized to A'(Civ). The data points are the AOD measured column 
densities of absorption redward of G2's systemic velocity, N(x) red , shown 
in Table [T] (column 8). We allowed the gas temperature and metallicity to 
vary for a range of hydrogen densities, n H . The data constrain the gas tem- 
perature to be 185,OOOK and limits log« H > -2.4. In Figure [TTl we provide 
the constraints on the warm phase metallicity. 



these ions are not part of the cool gas phase and are likely part of 
a separate "warm/hot" collisional ionized gas phase. In Figure [6] 
note that the Civ and Ovi have different kinematics that the lower 
ionization species further indicated that the absorption blueward of 
the galaxy systemic velocity is probing two gas phases. In the next 
section we model the warm gas. 



5.3 Warm Gas Phase 

We employed our own photo+collisional ionization code 
dChurchill & Klimekl 120121) to model the warm gas phase since 
it is optimally designed for optically thin gas with no ionization 
structure. In short, the code incorporates photoionization, Auger 
ionization, direct collisional ionization, excitation-autoionization, 
photo-recombination, high/low temperature dielectronic recom- 
bination, charge transfer ionization by H + , and charge transfer 
recombination by H° and He . All metal transitions and ionization 
stages for elements up to zinc are modeled . So lar abundance 
mass f ractions are obtained fromlDrainel d201 ll) and lAsplund et al.l 
(2009) and a lHaardt & Madaul ( 1201 lh ionizing spectrum is used 
for the ultraviolet background. The model inputs are the hydrogen 
number density (« H ), kinetic temperature and the gas metallicity 
while the model outputs are the electron density, and the ionization 
and recombination rate coefficients, ionization fractions and the 
number densities for all ionic species. Our models are consistent 
with those of Cloudy for log(n II )> -3.75. For log(n H )< -3.75, 
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Figure 11. — (left) The data points are the AOD measured column densities of absorption redward of G2 's systemic velocity, N(x)' €cl , shown in Table [JJ 
(column 8). — (middle) The photo+collisional ionization models obtained from Churchill & Klimek 12012) for Cm, Civ and Ovi using a fixed temperature 
of 185,000 K and [X/H] = -2.3. The spread in the models is due to the error in the measured Hi column density of log[N(Hi)]=16.03±0.18 (dotted lines). — 
(right) The allowed metallicity and n H for the model on the left. While n H is not well constrained, we find -2.50 < [X/H] < -1.93 and log;i H > -3.3 for the 
lower metallicity limit and logn H > -2.1 for the upper metallicity limit. 



cosmic ray heating in Cloudy tends to yield over-ionized clouds 
with large sizes (< 100 kpc). When cosmic ray heating is turned 
off, Cloudy has convergence issues and tends to produce ratios of 
column densities for adjacent species that do not follow a trend 
physical with ionization potential. We have reconciled these issues 
in our models and a deta iled comparison will be presented in 
IChurchill & Klimekl ( 120121) . 

5.3.1 Warm Gas Phase Blueward of Systemic Velocity 

In modeling the cold gas phase we noted that we are un- 
able to account for the measured logN(Civ)=14.11 and the 
logN(Ovi)=14.33. Thus, these ions are not part of the cool gas 
phase and are likely part of a separate warm collisionally ionized 
gas phase. Here we model the diffuse warm gas found blueward 
of G2's systemic velocity. Furthermore, all of the measured Cm is 
well modeled as being associated with the cool gas phase, however, 
there could possibly exist a small fraction of a measurable warm 
Cm component hidden. If we assume that the warm component is 
within the 0.2 dex measurement error of Cm, then there can be at 
most logA'(Cm)< 12.87 in the warm phase. 

In Figure|9]we show column density ratios, which is indepen- 
dent of Af(Hi), as a function of temperature and hydrogen density. 
The solid line is the Civ/Ovi column density ratio along with its 
measured error (dotted lines). The long and short dashed lines are 
the Cm/Civ and Cm/Ovi column density ratios, respectfully. The 
dashed curves are lower limits on the allowed temperature range 
for a given hydrogen density since logA'(Cm)< 12.87. Note that 
full self consistency is when the dashed curves fall below the solid 
curve. The models indicated that the gas has a hydrogen density of 
logn H > -3.5 and the Civ/Ovi provides temperature constraints of 
log? = 5.23 - 5.29, where the exact temperature is dependent on 

In Figure [9] we also show the gas metallicity as a function of 
A'(Hi) and n H . We are unable to measure the Hi column density 
associ ated with this warm component. However, Church ill et al.l 
d2007f) states that if an additional hydrogen component having 
logA'(Hi) > 16.8 was added to the total hydrogen column density, 



then it would measurably modify the shape of the Lyman limit. 
Thus, this provides a conservative estimate of the additional A'(Hi) 
that could be associated with the warm component since we assume 
it can not contribute to the Lyman limit. The column density limits 
provides a metallicity limit of [X/H] < -0.5. Note that the metal- 
licity scales with hydrogen column density by the same amount. 
It is likely that the metallicity is much lower, though we can not 
provide any further constraints. 



5.3.2 Warm Gas Phase Redward of Systemic Velocity 

Here we model the diffuse warm gas found redward of G2's sys- 
temic velocity. We allowed the gas temperature, the gas metallicity 
and the hydrogen density to vary to match the column density mea- 
surements, N(X) red , from Table [JJ (column 8). We only detect Cm, 
Civ and Ovi redward of the galaxy systemic velocity. The Hi as- 
sociated with the warm gas phase was computed from the AOD Hi 
column density redward of the galaxy systemic velocity. As shown 
in Figure [4] the Hi absorption profile is not saturated redward of 
the galaxy systemic velocity and contains a logN(Hi)=16.03±0.18. 

Figure [TO] shows the models column densities normalized to 
that of Civ. The data are well constrained with the warm phase 
gas having a temperature of T =185,000 K while there seems to 
be a larger range in the hydrogen density such that logn H > -2.4 
for a fixed metallicity. In order to accurately determine the range 
of metallicities and « H , we show the predicted column densities 
and metallicities as a function of n H in Figure [TT] Again, the mea- 
sured (left) and modeled (middle) column densities are shown. The 
model spread, indicated by the dotted lines, is dominated by the er- 
ror in the measured Hi column density. We find that the models do 
not well constrain the n H but provide reasonable constraints on the 
metallicity. We find -2.50 < [X/H] < -1.93 and log« H > -3.3 for 
the lower metallicity limit and logn H > -2. 1 for the upper metal- 
licity limit. 
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Table 2. Modeled Absorbing Gas Properties 



Name 


velocity [km s 1 ] 


logW(Hi) 


log(T) 


log(C/) 




l°g(« H ) 




[X/H] 






cold 


-240 < v abs < -3 


18.3±0.3 


3.82-5.23 


-3.25 <log({/) < 


-3.21 


-1.87<log(n H K- 


-1.83 


-1.68 « [X/H] 




-1.64 


warm blue 


-240 < v abs < -3 


<16.8 fl 


5.23-5.29 


log(C/) > 


-3.5 


log(n H ) > 


-3.5 


[X/H] 


< 


-0.5" 


warm red 


- 3 < Vabs < 220 


16.0±0.18 


5.27 


log(f/) > 


-3.3 


log(n H ) > 


-3.3 


-2.50 < [X/H] 




-1.93 



" See text in Section |5.3.1| for a discussion of the /V(Hi) upper limit used and how that affects the absorption metallicity. 



6 DISCUSSION 

The galaxy G2 has a stellar metallicity greater than solar, yet it con- 
tains [X/H] 2 halo gas detected in absorption at 104 kpc from 

the galaxy center. The data likely suggests that the absorbing gas is 
tracing cold accretion. A summary of the absorption gas properties 
is presented in Table [2] The kinematic difference^ seen between 
the high and low ionization states, as shown in Figure [6] and the 
photo+collisional ionization models suggest that the absorption is 
multi-phase. This may imply that the warm and cold gas physically 
arises in different locations or that the absorbing gas is not well 
mixed. 

The quasar sight-line probes absor ption within three de - 
grees of the projected minor axis of G2 jKacprzaketaT1l2011bh . 
which is a location rece ntly interpreted to be favorable for prob- 
ing galactic-scale win ds dBouche et alJl20lH ; lBordoloi et alj|201ll ; 
Kacp rzak et alj|2012h . G2 is dominated by a relatively old stellar 
population, yet it has measurable [On] emission (3.0 A) and thus is 
likely undergoing some current star-formation that could possibly 
produce winds. 

In Figure [7] we show a wind model that can account for the 
bulk of the warm gas kinematics, both redward and blueward of 
G2's systemic velocity, yet it does not reproduce the kinematics 
of the cold gas. This might imply that some of the warm gas 
could originate from winds, however, the 2 - 2.5 order of mag- 
nitude difference between the galaxy stellar metallicity and the 
absorption metallicity is inconsistent with a wind model (e.g., 
Oppenheim er et all |2010|) . Assuming a constant wind speed of 
150 km s , as suggested by the models, the absorbing gas would 
have left the galaxy 670 Myr ago. Given such a short time scale, 
G2 still expected to have solar metallicit y at that epoch acc ord- 
ing to mass-gas metallicity relations (e.g.. | Savaglio et al .120051) and 
simulations (e.g. JOppenheimer et alj2010[). The recyc ling for wind 
material is less than 1 Gyr (Oppenhei mer et ai] |2010h . thus, if the 
gas were wind material and was well mixed, it would be expected 
to have a metallicity similar to that of the parent stellar popula- 
tion. However, the efficiency at w hich the gas is mixe d within the 
halo environment is still unclear. dSiiacki et al.ll201 lh shows that 
poor gas mixing in simulations may be an artifact the suppres- 
sion of dynamical fluid instabilities demonstrated by comparing 
smooth particle hydrodynamics codes to that of the moving mesh 
code AREPO. However, some L LSs exhibit metallici ty variations 
across the absorption profile (e.g. JProchter et alj20 10). Thus, if the 
gas probed by this absorption system is poorly mixed, it is plausi- 
ble that the gas is entrain ed in a outflowing wind (see Alfiz as et al.l 
l2012tlSchave et alj|2007l) . Although, the current observational ev- 
idence suggests that Mgn entrained in outflowing material may, 



on average, have a m aximum projected extension of ~ 50 kpc. 
jBordoloi et al.l201lt) . Thus, the combination of the low metallicity, 
the relative kinematics between the galaxy and the absorbing gas, 
and the large projected distance between the gas and the galaxy, 
corroborate an accretion scenario. 

In summary, although the kinematic model works well for the 
warm gas, it is inconsistent with the expected outflow metallicities 
and this gas would have already been recycled several times and 
should have near solar metallicity. This leads to the idea that this 
meal poor gas could be infalling towards the galaxy. 

Galaxy G2 has a stellar mass of M, = 1x10" M Q and a viral 
mass of M vir = 8 x 10 12 M e . Cold mode accretion is less favored 
for galaxies with M yK > 10 12 M Q since the cold gas can become 
shock-heated as it enters the halo and ho t mode accretion becomes 
the dominant mode of accretion (e.g , Dekel & Birnboim 2006; 
Keres et alj|2009l; IStewart et alj|201 laM ; Ivan de Voort et aljboi ll; 



4 iDing. Charlton. & Churchill 1 2005) modeled the Civ and lower ioniza- 
tion states of Mgn, Mgi and Fen with the same kinematic structure. It is 
clear from the fit residuals that the kinematics are not well modeled with 
similar velocity structure although they do have some similarities. 



van de Voort & Schavell201lh . For massive galaxies, infalling gas 
can be shock heated to temperatures above 10 6 K, however, ac- 
cretion via dense filaments can be maintained past the shock re- 
gions, although have dramatically decreased covering fractions 
compared to low mass galaxies, and still accret onto the galaxy (e.g, 

i 1 — n i — 1 1 — — in i • — — \ 

Keres et al. 2009; Stewart et al. 201 la b; Faucher-Giguere et al. 
1201 ll ; Ivan de Voort et alj|201 ll ; Ivan de Voort & Schavell201 ll) . The 
hot gas pressure compresses cold streams and provides an efficient 
means of bringing col d pristine gas to the ho st galaxy and should 
be traced by LLS (e.g.. lFumagalli et alj|201 lh. A log MHi) = 18.3 
corresponds to Ap/p > 1000 at z < 1 (e.g. jDave et ai]| 19991) indi- 
cates that these high densities could make its way into the galaxy 
center. 

Although massive galaxy halos are built up from hot ac- 
cretion, comprising of 80 - 90% of the total accretion, it be- 
comes less important for accretion on to the ISM [50 - 70% 
hot m ode] dvan de Voort et al1l201 ll) . However. Ivan de Voort et ail 
d201 lh states that the cold gas accreting onto the ISM in massive 
galaxies increases with increasing simulation resolution. This could 
arise since gas clouds could be easily disrupted in SPH simula- 
tions, and the cold fraction would increase if higher densities were 
reached in higher resolution simulations. A lower hot fraction of 
gas accreted onto the galaxy occurs because the hot gas tempera- 
ture increases with the vira l temperature, or halo mass, resulting in 
longer cooling times (e.g.. lWiersma et al. 2009) and yielding less 
hot gas cooling to ISM temperatures. However these results could 
also b e dependent on the type of simulations used dSiiacki et al.l 
l201lh . 

Using the cosmolo gical formulae for virial quantities of 
iBrvan & Normanl dl998l) . we derived the radius, R vil = 380 kpc, 
the circular velocity v circ = 280 km s~' and the temperature, T vil = 
3 x 10 6 K. The quasar is probing along the minor axis of the galaxy 
within the viral radius at R/R vn ~ 0.3, well within the shock radius. 
Furthermore, the halo gas is expected to be heated to the viral tem- 
perature. The temperatures deduced from our models for the warm 
gas phase is 1.8 x 10 s K for gas blueward of the systemic velocity 
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and 1.9 xlO 5 K for the gas redward of the galaxy systemic velocity: 
both are an order of magnitude cooler than that of the viral tempera- 
ture. However, this gas can possibly cool after it had been shocked, 
although as mentioned above, the cooling time for hot shock-heated 
gas is much longer for massive galaxies. 

We can estimate a cool gas phase temperatu re by using the 
Dopple r parameters derived from Mgn absorption dKacprzak et al.l 
1201 lbi) . The Mgn Doppler parameters range from 2.1 - 10.8 km s -1 
and if we assume that all the broadening is thermal, we can compute 
upper limits on the temperature of 6.6x 10 3 - 1 .7x 10 5 K. These tem- 
peratures are upper limits since some of the line-broadening could 
be due to turbulence, and the Voigt profiles fits to the data assume 
a minimum number of "clouds"; if more clouds were inserted, than 
the velocity width of each line would decrease. Thus, the cold gas 
has temperatures well expected for cold mode accretion and not 
post-s hock heated cooling gas accre ting within a viral radius ratio 
of 0.3 l lvan de Voort & Schavell201 lh . 

The absorption kinematics also hints to an origin of cold- 
mode accretion. The metal poor cold-phase, and redward warm- 
phase, has kinematics consistent with extended disk rotation 
dKacprzak et al .l2010d ; ISteidel et alj2 002). This result has been in- 
terpreted using cosmological simulations to be a signature of cold 
gas accreting via filaments that drive the angular moment of the 
galax y, thereby mimicking i t s rotation out at large r impact param- 
eters dKacprzak et alj|2010al;lstewart et al1l201 lbl) . This is consis- 
tent with Ivan de Voort & Schavd d201 lh that shows cold accretion 
gas having a higher radial velocities and scales with increases mass 
compared to the flat radial velocity distribution of the hot-mode 
accretion as a function of mass. 

As we previously mentioned, the metallicities of the warm 
([XI H] = -2 to -2.5) and cold ([X/H] = -1.7) gas are low and 
appear to be consistent with metallicities expected for cold-mode 
accretion at R/R vlI ~ 0.3 dvan de Voort & Schave 201 1) and do not 
mimic the metallicities expected for th e ISM, the hot halo, and the 
host galaxy (Oppenheimer et al. 2010). The temperatures and kine- 
matics of the absorbing gas is also consistent with what is expected 
for cold-mode accretion. Thus, it is likely that the absorption is 
probing cold pristine gas infalling towards the center of the disk, 
further fueling star formation. 



7 CONCLUSIONS 

In this paper, we present detailed photo+collisional ionization mod- 
els and kinematics models of the multi-phase absorbing gas, de- 
tected within the HST/COS, HST/STIS, and Keck/HIRES spectra 
of the background quasar TON 153, associated with star-forming 
spiral galaxy at z = 0.6610. The sightline probes the projected mi- 
nor axis of the galaxy at projected distance of 0.3 virial radii, well 
inside the virial shock radius predicted for a galaxy of this mass, 
implying that if the gas is infalling that it is post shock heated ac- 
cretion or a cold filament. We obtained followup HST/COS data to 
study other metal-lines in order to determine the halo gas proper- 
ties and their origins. This galaxy was targeted as a candidate cold 
accretion p robe suppor t ed by k inematic and orientation results pre- 
sented bv lsteidel et al] d2002h ; lKacprzak et alj d2010all201 lbh . 
Our main results can be summarized as follows: 

(i) From g'r'i'Ks photometry and stellar population models, we 
determined that G2 is dominated by a ~ 4 Gyr stellar population 
with slightly greater than solar metallicity abundance and formed 
at redshift z ~ 2. We estimate an M, = 1 x 10" Mo implying an 
logM vir = 12.9. 



(ii) The low ionization states, Mgi, Sin, Mgn and Cm, have sim- 
ilar absorption kinematics, abundance ratios across the profile, and 
trace the bulk of the hydrogen, while Civ and O vi trace some of the 
same gas, their kinematics, abundance ratios, and their relative ab- 
sorption strengths differ. We infer that the low and high ionization 
states trace different gas phases. 

(iii) Modeling the cold gas blueward of G2's systemic velocity, 
N(X) bh,e , we constrain logr = 3.82-5.23,-3.25 «log(E/) < -3.21, 
-1.87 < log(n H )< -1.83, and -1.68 < [X/H] < -1.64. The gas is 
cold and very metal poor: consistent with cold accretion. We are 
unable to account for the measured N(Civ) and N(Ovi) when mod- 
eling the cold phase, thus, these ions are likely part of a separate 
warm collisionally ionized gas phase. 

(iv) A lagging halo model can account for all of low ionization 
absorption, hinting that this gas is coupled to the disk and simula- 
tions interpret this as a detection cold mode accretion. 

(v) Modeling the warm gas blueward of G2's systemic velocity, 
N(X) b ' ue , we find n H > -3.5 and logr = 5.23 - 5.29. Armed with 
only a conservative limit of hydrogen column density that could be 
associated with the warm component [A'(Hi) > 16.8], we estimate 
[X/H] < -0.5, although it is highly likely that the metallicity is 
much lower. 

(vi) Modeling the warm gas redward of G2's systemic veloc- 
ity, N(X) red , we find hot and metal poor gas with T =185,000 K, 
-2.50 < [X/H] < -1.93 andn n > -3.3. 

(vii) The quasar line-of-sight passes along G2's minor axis and 
a wind model can account for the observed Civ and Ovi redward 
and blueward of the galaxy systemic velocity. However, given the 
2 - 2.5 order of magnitude difference between the galaxy stellar 
metallicity and the absorption metallicity demonstrates the gas can 
not arise from galactic winds. 

It remains plausible that this low metallicity gas arises from 
unidentified satellites around the host galaxy or from the incom- 
plete mixing between metal enriched and metal poor halo gas, 
However, the combination of the relative kinematics, temperatures, 
and relative metallicities allows us to conclude that the multi-phase 
gas detected in absorption likely arises from cold accretion around 
this massive galaxy. For high mass galaxies the cold accretion 
cross-section is expected to be a few percent, so our absorption 
system and others cited in the literature could be a by-chance low 
probability inters ection of a filament, or the re solution effects in the 
simulations ("see lvan de Voort & S chave 2011) are underestimating 
the covering fraction of cold flows. This system also contradicts 
current results that predict that all absorption detected in quasars 
probing gas alon g the projected minor axis of galaxies is produced 
by wi nds dBordoloi et alj20ll] : lBouche et aljioT l ; Kacprza k et al.l 
l2012h : This is clearly not the case here. 
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